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SUBSTITUTE SPECIFICATION 

TITLE OF THE INVENTION 

COMPUTER-AIDED METHOD FOR PARALLEL CALCULATION OF THE OPERATING POINT 

OF ELECTRICAL CIRCUITS 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to German Application No. 199 
27 301.4 filed on June 15, 1999 in Gennany. and PCT Application No. PCT/DEOO/017754 filed 
on May 30, 2000, the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The computer-aided simulation of electrical circuits has attained increasing importance 
in the development of very large circuits, that is to say circuits having a very large number of 
elements. It is particulariy in the development of computer chips having a multiplicity, for 
example several hundred thousand transistors, that serial processing for the determination of 
the circuit quantities by a computer has been found to be unusable because of the excessive 
time consumption. 

[0003] In WO 98/24039, therefore, it is proposed to partition a large circuit and to have the 
partitions processed by different computers in each case. 

[0004] In the calculation, the operating point, that is to say the potentials of all nodes is usually 
determined first as basis for further analyses such as. for example, transient or alternating- 
current analyses. 

[0005] For the parallel calculation, an implementation of the Newton method is proposed in U. 
Wever. Q. Zheng et al.: "Domain Decomposition Methods for Circuit Simulation" (Proceedings 
of the 8th Workshop on Parallel and Distributed Simulation, PADS '94 Edinburgh, Scotland, UK, 
pp. 183-186, July 1994) and in U. Wever, Q. Zheng: "Parallel Transient Analysis for Circuit 
Simulation" (Proceedings of the 29th Annual Hawaii International Conference on System 
Sciences, pp. 442-447. 1996). The disadvantage is that convergence can only be achieved 
here when sufficiently good estimates of the operating point are available, due to poor 
convergence characteristics. As a rule, however, such good estimates can be achieved with 
difficulty or not at all in the case of large circuits. 

SUMMARY OF THE INVENTION 
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[0006] One aspect of the present invention is, therefore, based on the object of creating a 
computer-aided method for parallel calculation of the operating point of electrical circuits which 
ensures simple, reliable and fast calculation of the operating point of the circuit. 

[0007] In the method, the charging method, which has very good convergence characteristics, 
is used for the parallel calculation of the operating point. 

[0008] As described, for example, in H. Spiro: "Simulation of integrated circuits" (2nd edition, R. 
Oldenbourg Verlag Munich, Vienna 1990), the charging method uses the dynamic elements 
such as capacitances and/or inductances existing in an electrical circuit for calculating the 
operating point via a pseudo-transient analysis. In this process, the dynamic elements are 
charged up step-by-step to a value which, in principle, can have any magnitude, particularly to 
the value "1", starting from a value zero. 

[0009] In an embodiment according to the invention, dynamic elements such as, for example, 
capacitances or inductances having predetermined known values can be additionally installed at 
suitable places, for example at or between a number of nodes. This makes it possible to avoid 
problems in the calculation due to the number of dynamic elements being too small. 

[0010] For this purpose, the dynamic elements additionally inserted into the circuit can be 
reduced step-by-step to zero from an initial value, for example tending to infinity, for which the 
operating point calculation is trivial, so that the original circuit is again simulated. 

[001 1] In a preferred embodiment of the invention, at least one node, that is to say a junction of 
at least two current paths, of the circuit is connected by a capacitance to a predetermined value 
having a predetermined potential. 

[0012] In a further development of the invention, however, a capacitance, the second terminal 
of which is in each case connected to a predetermined potential, for example to ground, can 
also be connected to each node via all partitions. This procedure has the advantage that the 
calculation of the operating point, i.e. of the respective potentials of the individual nodes, is 
trivial for the circuit for an initial value for the capacitances which tend toward infinity, due to the 
equation, which is then explicit, being solved. Changing the value for the capacitances step-by- 
step then makes it possible to change the circuit simulation by suitable new selection of the 
value for the capacitances until a calculation of the operating point of the circuit is obtained for a 
value of the capacitances tending toward zero or almost toward zero. 
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[0013] By also implementing this charging method for parallel calculation of an electrical circuit, 
a very large circuit having a multiplicity of transistors can be advantageously calculated in a fast 
and simple manner even with a smaller number of dynamic elements. Various procedures are 
conceivable for the respective new determination of the value for the capacitance with the aim 
of allowing this to go toward zero, the decision criterion being the degree of difficulty of 
calculating the operating point of the respective preceding step. 

[0014] The operating point obtained by solving a nonlinear equation for a particular value of the 
capacitances in each case can be solved, for example, iteratively by the Newton method. For 
choosing the next value for "C", the number of iteration steps necessary for the solution can 
then be used until the value drops below a predetermined value for "C". 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] These and other objects and advantages of the present invention will become more 
apparent and more readily appreciated from the following description of the preferred 
embodiments, taken in conjunction with the accompanying drawings of which: 

[0016] In the text which follows, the invention will be explained by means of an exemplary 
embodiment shown in the drawing. In the drawing, the single figure shows a flowchart which 
represents the individual method steps of the method. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0017] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 

[0018] According to the drawing, a simulation of an electrical circuit is the starting point 1, for 
example in the circuit description language SPICE. This is partitioned in a first processing 
process 2, as described, for example, in WO 98/24039 as "Clustering method", so that individual 
partitions or parts of the circuits are obtained which can be calculated with the same degree of 
difficulty, if possible. 

[0019] Naturally, other applicable variants of partitioning are conceivable such as, for example, 
the "ratio cut method" described in N. Frohlich. B. Riess, U. Wever, Q. Zheng: "A new 
approach for parallel simulation of VLSI circuits on a transistor level" (IEEE Transaction on 
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Circuits and Systems - 1: Fundamental Theory and Applications, Vol. 45, No. 6, June 1998, 
pages 601 to 613), or even an arbitrary "manual" division. 

[0020] In a further process step 3, a grounded capacitance is added to each node of the circuit, 
that is to say to a junction of at least two conductors or current paths, respectively. Naturally, it 
is also conceivable to add a capacitance, the second terminal of which is connected to a 
predetermined potential, in which case both potentials and values for the capacitances of each 
node can differ. For reasons of efficiency, stability and economy of calculation, each 
capacitance is connected to the same potential with its second terminal, for example ground, an 
identically high value CO also being selected for all capacitances in step 4 for the same reasons. 

[0021] For this value CO, the operating point is then calculated for each partition or part-circuit 
in a further step 6, in which the required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being exchanged and included in the 
calculation of the operating points of adjacent partitions. 

[0022] In this process, a partition which will be called "master"' in the text which follows can 
advantageously take over the control of the charging process for reasons of efficiency. The 
master then determines the initial value CO for the capacitances, and it is also conceivable to 
predetermine the initial value externally, for example by the user. This value CO is then 
transferred to all other partitions, called slaves in the text which follows. Following this, the 
operating point is calculated both in the master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the energy sources existing in the 
circuit, is present at the circuit. 

[0023] Starting from an overall problem to be solved 

X = CiTJ, s,, . . . , Sp) 



where m is the unknowns of the master, 

p is the number of partitions. 

Si is the unknowns of partition i, 
the charging process leads to the system of differential equations, 

f(x,t)-fD^ = 0 
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where 

D = diag (C, ... C, 0, .... 0) and 
t is the time. 

[0024] This system of differential equations can be solved, for example, by the implicit Euler 
method which leads to the non-linear equation 

f(x''*',t''*') + iD(x''*'-x'') = 0 

where is the solution at time t'' for k = 0, 1 , etc and h is the step t''*^ - 1\ This non- 
linear equation can be solved, for example, iteratlvely by the Newton method 

xi::; =xr -(^f,(x-,t-)+^Dj' .[f(x-.t-)+^i>(xr -x^)] 

for n = 0, 1, etc. 

[0025] In the parallel calculation, the master then only calculates the unknowns "m", the 
unknowns "Sj" representing fixed areas for the master. To be able to calculate a solution of this, 
the slaves or slave partitions calculate the corrected values sf^^ in each iteration step of the 
master and report these to the master. For this purpose, the slaves must solve the system of 
non-linear equations. 

f,(sr',mr',t''-)+lD,((sr',mr)-(sNm'')) = 0 
n 

where fj and Dj are the corresponding equations and matrices for the slave "i". This non- 
linear equation can also be calculated by the Newton method and not all iteration steps have to 
be necessarily performed to convergence. 

[0026] The results of the slaves are then inserted into the system equations of the master 
whereupon the master can calculate the expression m^;^^ . 

[0027] In a next step 7, a new value "Cnew" for "C" is determined by the master and the slaves 
can also make suggestions. To achieve the aim, namely a value for "C" tending to zero or, 
respectively, less than a predetermined value e, and thus a calculation of the original circuit, the 
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choice of the new value "Cnew" for "C" is made dependent on the difficulty of calculating the 
preceding step with the value "Cow" for "C". 

[0028] In this context, various procedures are conceivable, for example the choice of "Cnew" in 
dependence on an analysis of the number of iteration steps which the master needed for 
calculating the non-linear system of equations for the preceding value "Cow" for "C": 

C / 

IS n < n, 



^new ~ 



Cold isn, <n<n2 
C , -2 is > n 

old ^ 



where "n" is the number of iterations of the master and "ni" "n/ are the parameters 
predetermined by the user. Naturally, instead of halving or doubling "Cow", respectively, other 
strategies for reducing or enlarging "Cow" are also possible. Additionally, the number of iteration 
steps needed by the slaves or slave partitions for solving their non-linear system of equations 
can also be taken into consideration, for example through the choice of 



n = max 



n, ^ iii , slavel,. . ^ , slavep 



V 



i=l 



where nj, siavej is the number of iteration steps of slave j. during the i-th iteration of the 
master. It must be noted in this context that the choice of a value for "C" only affects the 
efficiency of the method not the operating point itself. 

[0029] After a value for "C" of less than or equal to a predetermined value e has been reached, 
the calculation is terminated at branch 5 and the value for C can be set to "0" in a last step. The 
operating point of the original circuit is thus found as result 8 and can be output via output units 
such as for example, a screen, printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 

[0030] Due to the method described above, an operating point of a very large electrical circuit 
can be advantageously calculated in parallel by a multiplicity of computers or processors and 
the disadvantage of known parallel types of calculation, namely lack of convergence with 
unfavorable starting values can be avoided. Due to the parallel calculation, the expenditure for 
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calculation with regard to iteration steps and setting a value for "C" step by step, can be kept 
within tolerable limits due to the distribution to a number of processors. 

[0031] The operating.point thus calculated is then used as basis for further analyses, for 
example the alternating-current analysis of a circuit. 

[0032] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention. 



